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I  INTRODUCTION 


The  dispersive  nature  of  the  ionosphere  causes  a  systematic  time 
delay  of  radio  signals  that  propagate  through  it.  Depending  on  the 
precision  required,  this  time  delay  can  be  compensated  by  using  ionospheric 
models  (Klobuchar,  1975)  or  dual- frequency  transmissions  (see  for  example, 
Spilker,  1978).  Under  disturbed  conditions  there  is  a  random  "excess" 
delay  that  cannot  be  compensated  and  a  distortion  of  the  transmitted 
waveform.  Such  effects  are  a  concern  in  precision  navigation  systems 
such  as  the  Global  Positioning  System  (GPS). 

The  phenomenon  is  common  to  the  propagation  of  laser  pulses  in 
turbulent  media  (Liu  and  Yeh,  1980)  and  radiowaves  in  the  interplanetary 
medium  (Woo,  1975).  Yeh  and  Liu  (1977a, b;  1979)  and  Liu  and  Yeh  (1979) 
have  developed  a  detailed  theory  for  computing  the  excess  delay  and 
waveform  distortion.  There  are,  however,  few  data  sets  against  which 
the  various  theoretical  predictions  can  be  tested. 

In  this  report  we  have  developed  a  theoretical  formulation  of  the 
coherence  bandwidth  problem  (see  Section  II)  that  differs  from  Yeh  and 
Liu's  method  in  that  (1)  it  fully  accommodates  angle  effects  in  aniso¬ 
tropic  media,  (2)  it  gives  an  explicit  formula  for  the  single-point, 
two-frequency  coherence  function,  and  (3)  it  uses  asymptotic  rather  than 
Taylor  series  expansions,  thereby  preserving  the  dependence  on  the  power- 
law  index.  We  also  note  that  our  results  do  not  depend  on  the  precise 
value  of  the  outer-  or  inner-scale  cutoff  wavenumbers. 

To  validate  the  theoretical  predictions  (Section  III),  we  have 
used  data  from  the  Wideband  satellite,  which  transmits  phase-coherent 
signals  at  S-band,  I^-band,  UHF,  and  VHF.  The  UHF  signal  is  a  comb  of 
seven  equally  spaced  frequencies  (413.02  ±  nil. 47  MHz,  where  n  =  0,  1, 

2,  3)  that  was  uniquely  designed  for  coherence  bandwidth  measurements. 
Indeed,  the  name  "Wideband  satellite"  was  coined  because  of  this  capa¬ 
bility.  [See  Fremouw  et  al.  (1978)  for  a  detailed  description  of  the 
experiment. ] 
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The  theory  and  experimental  data  are  in  excellent  agreement.  As  in 
the  theory  of  the  time  structure  of  intensity  scintillation  under  strong 
scatter  conditions  (Rino  and  Owen,  1980),  a  single  parameter  characterizes 
the  combined  effects  of  changing  perturbation  strength  and  propagation 
distance.  The  effect  of  a  more  shallowly  sloped  spectral  density  function 
is  to  produce  more  frequency  decorrelation  for  a  fixed  perturbation  level. 
This  is  clearly  evident  in  the  data. 

The  comparisons  of  measured  and  predicted  frequency  coherence 
functions  validate  the  theory,  but  they  do  not  immediately  relate  to 
observable  system  effects  such  as  delay  jitter  or  pulse  distortion. 

Thus,  we  have  used  the  seven  coherent  signals  from  the  Wideband  satellite 
to  synthesize  disturbed  pulses  and  then  measure  the  actual  delay  jitter 
(Section  IV).  The  results  are  in  good  agreement  with  theoretical 
calculations  based  on  the  single-point,  two-frequency  coherence  function 
derived  in  Section  II. 

To  introduce  the  principal  quantity  of  interest,  we  recall  that  the 
ionosphere  is  a  linear  but  time-varying  transmission  medium.  Thus,  it 
can  be  characterized  by  a  time-varying  transfer  function  h(t;f)  (Bello, 
1963).  If  a  signal 

S (t)  =  Re[v(t)  exp  f 2«if c t} ]  (I- 1) 

is  transmitted,  the  received  signal,  SQ(t),  admits  a  similar  representation 
with  v(t)  replaced  by 

v  (t)  =1  v(f)h(t;f+f  )  exp  f2rtift}df  (1-2) 

o  I  c 

J~eo 

where  v(f)  is  the  Fourier  transform  of  v(t).  One  assumes  that  v(f) 

0  for  I f (  >  F/2 ,  and  that  F  «  fc> 

2 

An  incoherent  receiver  effectively  measures  (|vQ(t)|  >.  By 
direct  computation  from  Eq.  (1-2)  it  can  be  shown  that 


<|vo(t)|2>  =Jjv(f+6f/2)  V*  (f-6f/2)<h(t;fc+t-6f/2)  h*  (t;fc+f-6f /2» 

X  exp  {,2jti8ft}  dfif  df  .  (1-3) 

Thus,  the  entity  that  characterizes  the  randomly  dispersive  effects  of 

the  trans ionospheric  channel  is  the  single- point,  two- frequency  correlation 

function 

R(6f;f)  =  <h(t;f+6f/2)  h*(t ;f-5f /2)>  .  (1-4) 

If  coherent  detection  is  used,  a  similar  expression  for  the  power 
at  the  output  of  a  correlator  followed  by  a  low-pass  filter  can  be 
derived,  provided  that  h(t;f)  is  essentially  constant  over  intervals  of 
duration  1/F.  This  is  generally  a  very  good  assumption.  Thus,  the 
ionosphere- induced  delay  jitter  should  not  be  significantly  different 
for  coherent  or  noncoherent  detection. 


II  THE  SINGLE- POINT,  TWO- FREQUENCY  COHERENCE  FUNCTION 


Following  Fante  (1978),  we  first  consider  a  freely  propagating 
wavefield  u(px,z^;f^)  that  subtends  only  a  narrow  range  of  scattering 
angles.  If  we  take  z^  =  0  as  a  reference  plane  and  apply  the  Huygens- 
Fresnel  principle. 


2«iz 

a 


( 

( 


-i? 


2 


u(px  +  5,0;^)  df 


(II- 1) 


where  k  =  2jt/\  =  2rtf/c.  The  subscripts  l  and  l  denote  coordinates  along 
and  transverse  to  the  propagation  direction,  respectively.  Since  we  are 
only  interested  in  a  single-point  correlation  function,  there  is  no  need 
to  use  the  general  formulation  developed  in  Rino  and  Fremouw  (1977). 

Assuming  "frozen"  irregularity  structures  (the  Taylor  hypothesis), 
it  follows  that  h(t;f^)  =  u(vJLt,z^  ;f ^) .  Thus,  by  substituting  Eq.  (II-l) 
into  Eq.  (1-4),  we  have 


<h(t;f1)  h*  (t;f2)> 


klk2 

(2* z.)2 

i 


2  kl  .2  k2 
2*«  2\ 


Ru(Z?,f1,f2)  df  df' 


( 1 1-  2  ) 


where  R(Z^,f^,f2)  is  the  two-frequency  spatial  coherence  function  of 
u(p^,0;f).  We  have  assumed  that  u  is  spatially  homogeneous.  This  a 
assumption  does  not,  however,  require  homogeneous  phase  variations.  By 
changing  the  integration  variables  in  Eq.  (II-2)  to  f '  -  f  and 

X  =  (f  +  §/)/2,  we  can  isolate  and  evaluate  the  integral  over  X.  After 
a  series  of  algebraic  manipulations,  Eq.  (II-2)  can  then  be  put  in  the 
form 
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<h(t ;f x)  h*(t;f2)> 


1 

2rti 


k2)Zji 


kl'k2  1 

(kl  "  k2)2J 


Ru(zJ;flff2)  dzf 


(H-3) 


Equation  (II-3)  is  a  completely  general  expression  for  the  free-space 
propagation  of  the  two-frequency  coherence  function.  To  evaluate 

-4 

Ry(A§  ;f ^jf2) j  the  Rytov  approximation  can  be  applied,  as  was  done  by 
Ishimaru  (1972)  and  Fante  (1978).  It  is,  of  course,  simpler  to  use  an 
equivalent  phase-screen  whereby  all  diffraction  effects  develop  in  free 
space.  Because  the  results  based  on  the  equivalent  phase  screen  model 
seem  to  preserve  all  essential  aspects  of  the  scintillation  phenomenon, 
we  have  used  that  model  herein. 


To  apply  the  phase  screen  model  we  let 


udiOjfj)  =  exp  fi^XAN^s)} 


(11-4) 


where  AN (£)  is  the  perturbation  to  the  integrated  electron  content.  If 
— *  ™ 

AN^(?)  is  gaussian,  it  is  readily  shown  from  Eq.  (II-4)  that 


( 

Ru<^  ;fpf2)  =  exp  I  "r 


1  -  e 


2  AN 


w 


l 


2,  2 
reX 


2e 


1  -  e 


(AN' 


;>! 


(II-5) 


where 


f  =  c/X  -  (fj  +  f2)/2 
6  f  =  f 2  -  fj 

e  =  6  f / (2f) 

and  (A§ )  =  ((AN^(f)  -  /^N^(5/))^)  is  the  structure  function  for 
AN^(E).  The  gaussian  assumption  should  be  regarded  as  a  sufficient 
not  strictly  necessary  condition  (Rino,  1979b). 

If  e  «  1,  then  Eq.  (II-5)  simplifies  to 


(II-6a) 

(II-6b) 

(II-6c) 


but 


Ru(^;fj8f)  =  exp  {-D^(A£)}  exp  f-2e2CT^} 


(II- 7) 
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2  -» 
where  0-  is  the  phase  variance,  and  D  (Z£)  is  the  phase  structure 
60  60 

function.  Hereafter,  it  is  understood  that  any  quantity  that  depends 
implicitly  on  frequency  is  to  be  evaluated  at  the  mean  frequency, 

Eq.  (II-6a).  Note,  however,  that  the  mean  frequency  f  need  not  equal  fc. 

By  similarly  approximating  (kj^)  /  (kj-k£) ,  changing  variables  in 
Eq.  (II-3),  and  substituting  from  Eq.  (II-7)  we  obtain  the  general 
result 


R(6f;f)  =^||exp  f-i  Z^2/2]  exp  f-Dg^(Z^[2Z  6f/f]1/2)}  dz£ 


X  exp  {-a5202e2} 


Since  we  are  working  in  a  transverse  coordinate  system,  D  (£%)  is 

60 

functionally  dependent  on  the  quadratic  form 


(II-8) 


2 

y 


C^2  -  B^x^y  +  A^2 
A7C7  -  B'2/4 


(II-9) 


where  A7,  B7,  and  C7  are  defined  by  Eqs.  (26a),  (26b),  and  (26c)  in 
Rino  (1979a). 

By  first  rotating  variables  to  remove  the  term  and  then 

performing  a  series  of  variable  changes,  it  is  possible  to  reduce  Eq. 
(II-8)  to  the  single  integral 


R(6f;f) 


=  i\/f 


8  - a 


CO 


(ury)  exp  {-icuB} 


X  exp  f-D  ([u)2Z  6f/f]1/2)}  duj  exp  f-a. .2e2} 
0  (5  50 

where 

A"  =  -j(A7  +  C7  +  D7) 

C"  =  -|(A7  +  C7  -  D7) 


(II- 10) 

(Il-lla) 

(II-llb) 
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(II- 11c) 
(II- lid) 


D'  =  \/(A'  -  c')2  +  B/2 
a  =  (A."  -  C")/4 
g  =  (A"  +  C")/4  .  (II- lie) 

The  same  transformation  was  used  in  Rino  (1979a)  [cf.  Eqs.  (29a),  (29b), 

and  (29c)].  For  isotropic  irregularities.  A"  =  c"  =  1,  so  that  p  =  0.5 

2 

and  <y  =  0.  For  highly  anisotropic  irregularities,  a"  ~  a  ,  where  a  is 

2  2 

the  axial  ratio,  and  C"  =  1.  For  large  a,  a  ~  (3,  although  p  -  o'  = 

A"  C"  /4  ~  a2/4. 

In  a  power- law  environment  the  general  form  of  the  phase  structure 
function  is  too  cumbersome  to  be  of  practical  value.  Thus,  approximations 
must  be  used  that  take  advantage  of  the  fact  that  the  outer  scale  is 
typically  much  larger  than  the  scale  sizes  that  make  a  significant  contri¬ 
bution  to  the  signal  over  time  intervals  of  interest.  Two  different 
approximations  have  been  used. 

For  a  three-dimensional  spectral  density  function  of  theform 

(q)  =  Csq"(2V+1>  (II'12) 

e 


with  0.5  <  v  <  1.5,  it  can  be  shown  by  asymptotic  methods  that 


Vy)  ~ 


|y|2v’1 


where  C  is  the  phase  structure  constant, 


i  2r(1.5-v) 

2n  r(v+0.5)(2v-l)22v_1 


(11-13) 


(11-14) 


In  this  model,  which  has  been  used  extensively  in  neutral  turbulence 
studies,  the  spectral  density  function  and  the  structure  function  have 
complementary  power-law  forms.  There  is,  moreover,  no  dependence  on 
the  inner-  or  outer-scale  cutoff  wavenumbers. 
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Alternatively,  one  may  use  a  Taylor  series  expansion  of  D  (y)  and 

60 

retain  only  the  quadratic  term.  In  that  case, 


°6/y) 


(11-15) 


where  the  form  of  is  given  in  the  Appendix.  Comparisons  of  the 
quadratic  and  asymptotic  approximations  are  also  given  in  the  Appendix. 
The  quadratic  approximation  is  strictly  valid  only  for  steeply  sloped 

spectral  density  functions  such  that  v  >  1.5.  It  should  be  noted  that 

2 

both  a  and  D  depend  on  the  outer-scale  wavenumber  q  ,  and  D  depends 
60  i  o  i 

on  the  inner-scale  cutoff  wavenumber  as  well. 

If  Eq.  (11-13)  is  substituted  into  Eq.  (II- 10)  the  result  is 


R(8i 


f;f)  =  2  \/@2  -  a2  I  J,(«w)  exp  f-icup } 

Jo 


X  exp 


u 


f 


v-0.5 


v-0.5)  .  ,2.2, 

co  dm  exp  2e  } 


(11-16) 


where 


H  =  G 


\2Z\ 


v-0.5 


(11-17) 


The  factor  G  accounts  for  anisotropic  media.  It  is  defined  and  discussed 

in  detail  in  Rino  (1979a, b).  The  quadratic  approximation  gives  rise  to 

a  similar  form  with  the  exception  that  v-0.5  is  replaced  by  unity  in  Eqs. 

(11-16)  and  (11-17),  and  C  \  in  Eq.  (11-17)  is  replaced  by  2<j,^D  q^. 

60  50  I  ° 

Unfortunately,  Eq.  (11-16)  cannot  be  analytically  evaluated  for  v 

within  the  admissible  range  0.5  >  v  >  1,5.  However,  following  Fante 

(1978),  who  encountered  a  similar  integral  in  this  context,  one  can 

replace  the  v-0.5  term  in  Eqs.  (11-16)  and  (11-17)  by  unity.  The  integral 

can  then  be  evaluated  analytically,  giving  the  result 
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R(6f;f)  a- 


2  r  2„  2, 

exp  f-p^e  } 


16  -  iH 


6f 

f 


)  -„2 


(11-18) 


We  shall  see  that  Eq.  (11-18)  predicts  a  correlation  level  that  is  too 
high  for  a  given  value  of  H.  Under  the  quadratic  approximation,  however, 
Eq.  (11-18)  is  exact. 

In  effect,  Eq.  (11-18)  becomes  a  better  and  better  approximation  as 
the  spectral  index  parameter  approaches  v  =  1.5.  For  more  steeply  sloped 
spectra,  one  expects  the  quadratic  approximation  to  give  accurate  results. 
The  crucial  factor,  therefore,  is  the  spectral  index.  We  shall  see  in 
Section  III  that  the  data  clearly  fall  below  the  predictions  of  Eq. 
(11-18),  thereby  further  supporting  the  more  shallowly  sloped  spectral 
density  functions  that  have  now  been  verified  by  a  number  of  different 
studies  (Rino,  1979a, b;  Rino  and  Owen,  1980). 


I 

l 

\  i 


4 


| 

i 

i 

III  DATA  ANALYSIS  I 

j 

As  noted  in  the  Introduction,  the  Wideband  satellite  transmits  seven 
equispaced,  phase- coherent  signals  at  UHF.  Thus,  21  values  of  (u(f i)u*(f^)) 
can  be  measured  with  different  values  of  f^  and  f It  is  convenient, 
however,  to  use  only  the  three  symmetric  frequency  pairs  about  413.0244 
MHz.  In  that  case  the  mean  frequency  f  =  (f^  +  f.^) /2  is  invariant  and 
identically  equal  to  the  center  frequency.  The  corresponding  values  of 
Sf/f  are  0.056,  0.111,  and  0.167. 

To  obtain  a  data  base  of  highly  disturbed  passes  we  have  used  the 
same  equatorial  Wideband  passes  as  were  processed  and  analyzed  in  Rino 
and  Owen  (1980).  The  seven  UHF  channels  were  detrended  with  a  10-s 
detrend  filter  cutoff  at  a  100-Hz  sample  rate.  For  single- time- point 
measurements  the  sample  rate  is  not  critical.  Thus,  to  save  computer 
time  the  lower  data  rate  was  used.  The  28  complex  correlation  coeffi¬ 
cients  were  then  computed  and  recorded  for  further  processing. 

Now,  from  Eqs .  (11-16)  and  (11-17)  we  see  that  R(6f;f)  depends  on 
the  propagation  geometry  through  the  a  and  B  terms  as  well  as  H.  However, 

for  the  highly  elongated  irregularities  near  the  geomagnetic  equator  we  i 

2  2  2  2  2 
can  use  the  limiting  values  a  ~  (a  +l)/4,  B  ~  (a  -l)/4,  and  p  -  a  ~  a  /4, 

where  a  is  the  axial  ratio.  Thus,  the  principal  variable  in  Eq.  (11-16)  ; 

is  the  combined  perturbation  strength  and  propagation  distance  parameter, 

H*  j 

We  also  note  that  Eq.  (11-16)  admits  a  weak  dependence  on  the  rms 

phase,  <j  ,  Since  the  value  of  the  rms  phase  itself  depends  on  the  i 

off  1 

interval  over  which  it  is  measured,  R(6f;f)  will  exhibit  a  similar  non¬ 
stationary  behavior.  From  Eq.  (II-5)  it  can  be  seen  that  the  phase  term 
is  a  remnant  of  the  purely  dispersive  behavior  of  the  ionosphere.  Its 
contribution  is  generally  small,  and  in  any  case  is  readily  computed 

because  we  measure  <r„  . .  Alternatively,  one  can  remove  linear  phase  trends 

off 

across  the  UHF  comb  before  the  correlation  coefficients  are  computed. 
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The  parameter  H  depends  on  the  geometrical  enhancement  factor,  G, 

the  phase  structure  constant,  C  ,  and  the  Fresnel  parameter,  Z.  Since 

2  2  1 

the  phase  turbulent  strength,  T,  can  be  measured,  and  T  «  GC  v  v“ 

o  50  Gif 

[see  Rino  (1977),  Eq.  (18)3,  we  could  determine  GC  by  first  estimating 

00 

the  effective  scan  velocity  v^^.  Hie  problem  is  that  v  depends  on 
both  the  assumed  layer  height  and  the  propagation  geometry  as  does  Z. 

If  this  approach  were  pursued,  therefore,  the  data  sorting  would  depend 
on  an  unknown  parameter — namely,  the  height  of  the  equivalent  phase  screen. 

We  note,  however,  that  under  conditions  of  weak  scatter  for  S^, 


s,  =  c  z 

4  p 


.v-l/zf  r  (2 ,5-v) /2 

r  (v+0.5)/2  (v-0.5)_ 


(III-l) 


where  J  ~  r(v)/[yit  F(v+l/2)]  for  large  a  (Rino,  1979a).  Since  G~  1  under 

2 

the  same  conditions,  H  and  S^  are  simply  proportional.  The  proportionality 
constant  depends  on  v  but  is  readily  calculated  by  comparing  Eqs.  (11-14), 
(11-17),  and  (III-l).  In  our  data  analysis  we  used  the  measured  L-band 
S^  scintillation  index,  which  was  generally  in  the  weak  scatter  regime 
(S^  <0.6),  to  compute  H  for  an  assumed  value  of  v. 

In  Figure  1,  measured  values  of  R(fif;f)  plotted  against  H  are  shown 
for  a  set  of  disturbed  passes  recorded  at  Kwajalein.  The  numbers  above 
each  data  point  indicate  the  number  of  measured  R(fit;t)  values  that  fell 
in  the  particular  H  bin.  The  error  bars  indicate  the  maximum  excursion 
for  the  measured  standard  deviation  of  the  Af/f  =  0.056  data  points,  and 
the  minimum  excursion  for  the  corresponding  deviation  of  the  Af/f  =  0.167 
data  points.  In  Figure  2  similar  measurements  for  the  Ancon  data  are 
shown. 

Both  data  sets  show  that  measured  values  of  R(Af;f)  are  well  ordered 
by  the  H  parameter.  The  exact  value  of  H  depends  on  the  v  value  and 
weakly  on  the  propagation  geometry.  The  particular  value  v  =  1.5,  which 
corresponds  to  a  one- dimensional  phase  spectral  index  of  2.5,  is  a 
representative  median  value  from  spectral  analysis  of  the  corresponding 
phase  data  (Rino,  1979a;  Rino  and  Owen,  1980). 
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FIGURE  1  MEASURED  SINGLE-POINT  TWO-FREQUENCY  CORRELATION  FUNCTIONS  FROM  KWAJALEIN  DATA  PLOTTED 
AGAINST  H  AS  DERIVED  FROM  SIMULTANEOUS  L-BAND  Sa  DATA 


FIGURE  2  MEASURED  SINGLE-POINT  TWO-FREQUENCY  CORRELATION  FUNCTION  FROM  ANCON  DATA  PLOTTED  AGAINST 
H  DERIVED  FROM  SIMULTANEOUS  L-BAND  S,  DATA 


To  compare  the  measured  R(Sf;f)  values  with  the  theory,  we  have 

2  2 

numerically  integrated  Eq.  (11-16)  with  a  =  (a  +  l)/4,  and  6  =  (a  -  l)/4, 
and  different  values  of  v.  Only  Af/f  =  0.167  was  used  because  it  gives 
the  maximum  decorrelation  for  a  given  value  of  H,  and  the  integral  is 
most  easily  evaluated  in  such  regimes.  In  Figure  3  the  theoretical 
results  are  plotted  together  with  the  approximate  formula  Eq.  (11-18). 

Also  shown  in  the  figure  are  measured  values  of  R(6f;f)  from  the  Ancon 
data. 

We  first  note  that  the  approximate  formula,  Eq.  (11-18),  always 
predicts  a  smaller  frequency  decorrelation  than  the  exact  formula,  and 
that  the  amount  of  decorrelation  for  a  given  value  of  H  increases  with 
decreasing  v.  The  actual  data  values  systematically  fall  below  the 
theoretical  curves  for  decreasing  v  values  as  H  increases.  This  behavior 
suggests  a  variable  spectral  index,  which  has  indeed  been  observed  and 
used  to  improve  the  theoretical  fit  to  intensity  coherence  measurements 
under  strong  scatter  conditions  (Rino  and  Owen,  1980).  The  spectral 
index  decreases  with  increasing  perturbation  strength.  Thus  the  theory 
and  data  are  in  excellent  agreement. 
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PHASE  SCREEN  THEORY 


IV  EXCESS  DELAY  JITTER  AND  PULSE  DISPERSION 


We  have  developed  a  model  that  characterizes  the  single- point  two- 
frequency  coherence  function  and  showed  that  it  is  in  good  agreement 
with  data  from  the  Wideband  satellite.  As  yet,  however,  we  have  not 
developed  a  simple  means  of  assessing  the  impact  of  the  randomly 
dispersive  effects  on  broadband  waveforms.  To  this  end,  Yeh  and  Liu 
(1977a)  have  computed  the  temporal  moments 


00 

M  =|  t  <|v 
n  I  n  '  o 

J- 


<t)|2>  dt  . 


(IV-1) 


They  showed,  for  example,  that  if  j’tv(t)dt  =  0,  then 
F/2 


T  A  M  = 

Td  =  M1  2jti 


Til  lJ<» 
J-  F/2 


iff  R(fif;?) 


df 


(IV- 2) 


J6f=0 


and 


F/2 

n  ■  Mo  -  I  t2lv(t)|2dt  =  --1-  -2  f 

1  (2*i)  J- F/2 


J1 


|v(f) 


Laef 


R(6f 


;f)  df. 

Jsf=o 


(IV-3) 


Mj  is  a  measure  of  the  average  delay  jitter  and  Eq.  (IV-3)  measures  the 

waveform  spreading.  However,  from  Eq.  (11-10),  the  nth  derivative  of 

R(6f;f)  is  proportional  to  1 6 f /f  |v”  (^^.S),  singular  at  fif  *  0 

when  0.5  «r-  v  <  1.5.  This  implies  that  in  a  power-law  environment 
2 

<[vo(t) I  )  develops  a  slowly  decaying  tail. 

The  effect  is  not  an  artifact  of  the  model,  but  rather  an  indication 
that  one  cannot  meaningfully  evaluate  Eq.  (IV-1)  over  an  arbitrarily 
large  time  interval.  The  singularity  can,  of  course,  be  removed  by 
introducing  an  inner-scale  cutoff.  We  believe,  however,  that  since  an 


inner-scale  cutoff  has  not  been  detected  in  scintillation  data,  it  is 
more  realistic  to  use  Eqs.  (IF-2)  and  (IV-3)  with  the  derivatives 
evaluated  at  §f  =  tc*  instead  of  §f  =0  is  on  the  order  of  the  waveform 
duration. 

It  happens,  however,  that  if  the  approximate  form  of  Eq.  (11-18)  is 
used,  the  singularity  is  also  removed.  It  is  then  readily  shown  that  if 
F  «  fc, 

Td~  -dr  2'P .  2  (IV-4) 

c  P  -  O' 

and 


(IV- 5) 


For  isotropic  irregularities  at  normal  incidence,  g  =  1  and  <y  =  0.  In 

2 

that  case,  -  (t^)  .  The  same  relationship  holds  in  a  deterministically 
dispersive  uniform  ionosphere  where  oc  N^/f^,  and  is  the  integrated 
content. 


It  remains  to  show  how  good  a  measure  Eq.  (IV-4)  or  more  refined 
estimates  are,  since  they  have  practical  ramifications.  We  first  note 
from  Eq.  (1-2)  that  the  modulation  imparted  on  a  sinusoidal  transmission 
is  h(t;f).  Thus,  from  the  Wideband  satellite  UHF  comb  of  seven  frequencies 
we  can  synthesize  a  pulse  as 


3 

v°(t)  =  ^  ^  h(t;fk)  exp  f2itikAfT]  (IV-6) 

k=-3 


where  =  413.02  ±  kAf,  and  Af  =  11.47  MHz.  For  example,  if  there  is  no 
disturbance,  so  that  h(t;f^)  s  1,  then  it  is  easily  shown  from  Eq.  (IV-6) 
that 

I* „<T>|2  .hsLumiI  .  (iv- 7, 

sin  [nAff] 
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The  pulse  is,  of  course,  non-causal  and  it  has  range  ambiguities 

at  multiples  of  1/Af  =  0.0872  ps.  This  is  of  no  consequence  for  our 

purposes  here,  however,  and  we  can  use  Eq.  (IV-6)  at  different  t  values 

2 

to  first  estimate  <|vq(t)|  )  and  then  compute  various  signal  moments  as 
defined  by  Eq.  (IV-1).  Because  the  procedure  is  somewhat  time  consuming, 
only  the  first  moment  was  actually  computed. 

2 

A  typical  set  of  (|vo(t)|  )  estimates  is  shown  inFigure  4.  The 
pulses  were  averaged  over  a  20- s  data  interval.  The  numbers  on  the 
right-hand  edge  of  each  pulse  estimate  give  the  corresponding  H  value. 

For  H  ^  3,  the  pulses  show  severe  distortion  as  we  should  expect  from 
the  data  in  Figures  1  and  2. 

Estimates  of  the  first  moments  together  with  the  theoretical  pre¬ 
diction  of  from  Eq.  (IV-4)  are  shown  in  Figure  5.  There  is  a  large 
amount  of  scatter  in  the  data  because  of  the  sensitivity  of  simple  moment 
estimates  to  small  fluctuations  at  large  t  values,  which  is  the  nature 
of  ionospheric  coherence  bandwidth  effects.  In  any  case,  the  crude 
theoretical  calculations  upon  which  Eq.  (IV-4)  is  based  seem  to  give  an 
acceptable  result  for  engineering  purposes.  There  is  a  tendency  for  the 
actual  delay  jitter  to  be  larger  than  the  predicted  value  from  Eq.  (IV-4). 
This  is  expected  since  the  observed  decorrelation  is  always  somewhat 
greater  than  that  predicted  by  Eq.  (11-18)  from  which  Eq.  (IV-4)  was 
derived . 
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FIGURE  4  SYNTHESIZED  PULSES  FROM  WIDEBAND  SATELLITE  DATA  SHOWING  EFFECTS  OF  COHERENT  BANDWIDTH  LOSS 


FIGURE  5  MEASUREMENTS  OF  DELAY  JITTER  FROM  SYNTHESIZED  PULSES  USING  WIDEBAND  SATELLITE  DATA  TOGETHER 
WITH  THEORETICAL  PREDICTION  FROM  EQ.  (IV-4) 


V  DISCUSSION 


In  this  report  we  developed  a  theoretical  model  for  computing  the 
single-point,  two- frequency  coherence  function  that  characterizes  the 
ionospheric  coherence  bandwidth  loss  (Section  II).  We  followed  the 
method  of  Fante  (1978),  who  proceeded  from  the  Huygens- Fresnel  principle 
and  then  used  the  Rytov  approximation  to  determine  the  input  conditions. 
In  our  approach  we  simplified  the  computation  by  neglecting  the  dif¬ 
fraction  effects  as  predicted  by  the  Rytov  method,  thereby  assigning 
all  the  diffraction  effects  to  free-space  propagation. 

As  discussed  in  detail  in  Rino  (1979a, b),  this  approach  generally 
gives  simple  equivalent  or  "lumped"  parameters  that  characterize  the 
ionospheric  perturbation  and  the  resultant  diffraction  effects.  Indeed, 
in  the  case  of  the  shallowly  sloped  spectra  (v  <  1.5)  that  seem  to 
characterize  ionospheric  irregularities,  a  single  parameter  H  defined  by 
Eq.  (11-17)  characterizes  the  combined  effects  of  increased  perturbation 
strength  and  propagation  distance.  The  H  parameter,  like  the  U  parameter 
that  characterizes  the  corresponding  effects  in  intensity  scintillation 
data  (Rino,  1979b),  does  not  depend  on  the  inner-  or  outer-scale  cutoff 
wavenumbers . 

The  theory  is  in  good  agreement  with  measured  frequency  coherence 
functions  using  the  Wideband  satellite  data  (Section  III).  A  crude 
approximation,  Eq.  (11-18),  to  the  integral  form  of  the  theoretical  model 
tends  to  underestimate  the  observed  amount  of  frequency  decorrelation, 
but  gives  acceptable  results  for  engineering  applications.  In  this 
regard,  we  also  synthesized  a  pulse  using  the  Wideband  comb  of  seven 
equally  spaced  UHF  transmissions  and  then  measured  the  delay  jitter. 
Moment  estimates  by  their  very  nature  produce  large  amounts  of  scatter, 
but  good  overall  agreement  between  the  measurements  and  the  theory  was 
obtained. 
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In  a  series  of  papers,  Yeh  and  Liu  have  developed  a  detailed 
formalism  for  computing  the  signal  moments  that  characterize  the  iono¬ 
sphere-induced  delay  jitter  and  pulse  smearing.  Their  method  relies  on 
a  Taylor  series  expansion  that  suppresses  the  dependence  of  coherence 
bandwidth  effects  on  the  power- law  index  and  depends  on  both  the  inner- 
and  outer-scale  cutoff  wavenumbers.  If  the  same  (quadratic)  approxi¬ 
mation  is  used  in  our  own  model,  Yeh  and  Liu's  results  can  be  recovered. 
For  example,  our  formula  for  Td  is  equivalent  to  t^  given  by  Eq.  (5)  in 
Yeh  and  Liu  (1977). 

To  show  the  differences  in  the  predicted  values,  in  Figure  6  we  have 
computed  using  both  the  small  qQ  (asymptotic)  and  quadratic  approxi¬ 
mations  for  respective  spectral  indices  in  ranges  where  the  approximations 
are  valid.  For  the  quadratic  approximation,  an  outer  scale  of  10  km  was 
used  (qQ/2jt  =  10  km),  and  an  inner  scale  of  100  m.  The  results  were 

computed  for  a  1-GHz  signal  and  plotted  against  C  to  facilitate  compari- 

s 

sons  with  Figure  8  in  Rino  (1979a). 

It  can  be  seen  that  the  quadratic  theory  predicts  much  smaller  delay 
jitter  than  does  the  small  q^  asymptotic  theory.  These  results  are  not 
in  conflict,  but  illustrate  the  critical  dependence  of  the  theory  on  the 
power- law  index.  All  the  Wideband  data  that  have  been  analyzed  to  date 
have  favored  the  more  shallowly  sloped  spectra  where  the  asymptotic  theory 
is  valid.  We  have,  moreover,  recently  found  evidence  that  the  spectral 
index  varies  with  changing  perturbation  strength  (Livingston  and  Rino, 
1980).  The  effects  of  a  varying  spectral  index  are  clearly  evident  in 
the  frequency  coherence  data  presented  in  Section  III. 

The  theory  developed  in  this  report  can  be  applied  essentially  over 
the  full  range  of  spectral  indices  that  might  occur  in  a  disturbed 
nuclear  environment.  The  simple  formulas  for  the  delay  jitter  and  pulse 
broadening  can  be  applied  to  get  a  rough  estimate  of  the  effects  in 
precision  navigation  systems,  for  example.  More  refined  calculations 
and/or  simulations  can  be  performed  by  using  the  full  frequency  coherence 
function. 
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Appendix 


APPROXIMATIONS  TO  THE  PHASE  STRUCTURE  FUNCTION 


The  general  form  of  the  phase  structure  function  corresponding  to 
the  general  power- law  model  summarized  in  Table  2  of  Rino  and  Fremouw 
(1977)  is 
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r  r(v-l/2) 

p  2rtT(v+l/2) 
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v-1/2 


(qQy/2) 


K 


v-1/2 


(2sf? T 


(qQy/2)^ 


(A-2) 


r  =  9  /q,  «  1 
o  1 

and 

*  ’  1  rV'1/2Kv-l/2<2r> 
-  r(v-i/2)  . 

A  plot  of  £(y)  is  shown  inFigure  (A-l). 

It  was  shown  in  Rino  (1979b)  that 


(A-3) 


(A-4) 


£<y)  ~  | y | 2v” 1  (a-5) 

as  long  as  0.5  <  v  <  1.5.  The  formal  Taylor  series  expansion  of  Eq. 
(A-2)  can  be  written 

oe 

£<y>  =  VV)2n  (A“6) 

k=l 


29 


FIGURE  A-1  PLOT  OF  NORMALIZED  STRUCTURE  FUNCTION  SHOWING  DEPENDENCE  ON 
SPECTRAL  SLOPE 


where 


D  =  * - "=- 

n  (2n) I  9y  n  y=0 


(A- 7) 


For  our  purposes  here  we  need  only  consider  the  first  nonzero  coefficient 


D 


1 


r-3/2  K 


v-3/2 


2  l°g  <2r> 


r(v-i.5) 

4r(v-0.5) 


<2e)/N 
v  =  1/5 

v  >  1.5 


(A- 8) 


Higher-order  coefficients  become  increasingly  more  sensitive  to  the 
inner-scale  cutoff. 

In  Figure  (A-2)  the  small  qQ  approximation  is  shown  on  an  expanded 
plot  of  the  exact  form  of  £(y).  In  Figure  (A-3)  the  corresponding  plot 
for  the  quadratic  approximation  is  shown.  The  asymptotic  approximation 
is  applicable  when  v  <  1.4  and  improves  as  v  decreases,  whereas  the 
quadratic  approximation  applies  when  v  >  1.5  and  improves  with  increasing 
v.  Both  approximations  are,  of  course,  only  valid  for  qQy  «  (2jt) 

Because  of  the  very  large  outer-scale  wavelength  that  characterizes 
ionospheric  irregularities,  however,  this  is  generally  a  very  good 
approximation. 
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